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Abstract 

Transition frequencies of two spherical gas bubbles interacting in an acoustic field 
are discussed theoretically. In the present study, transition frequency is defined as 
the frequency of external sound for which the phase difference between a bubble's 
pulsation and the external sound is tt/2. It is shown by a linear theory that a bubble 
interacting with a neighboring bubble has three (or fewer) transition frequencies but 
only two natural frequencies. This result means that the bubble has a characteristic 
frequency besides the natural frequencies. 

Key words: Two-bubble dynamics, Radiative interaction, Natural frequency, 
Phase reversal 

PACS: 43.20.+g, 47.55.Bx, 47.55.Dz 



It is known that radiative interaction between gas bubbles varies their natu- 
ral (or resonance) frequencies, around which the bubbles indicate resonance 
response. This variation in the natural frequencies has been the subject of 
many studies [1-13]. In where two spherical gas bubbles of different 

radii interact acoustically, a bubble has two (or one) natural frequencies (see, 
e.g., [2,5,6,8]). 

The aim of this Letter is to show that a gas bubble interacting with a neigh- 
boring bubble has an alternative characteristic frequency; we now call this a 
transition frequency. In the present work, a transition frequency is defined as 
the frequency of an external sound (the driving frequency) for which the phase 
difference between a bubble's pulsation and the external sound becomes tt/2. 
It is well known that, in a single-bubble case, the phase difference between 
the bubble and an external sound becomes tt/2 when the driving frequency is 
equal to the bubble's natural frequency [14]. However, it is shown in this Let- 
ter that, in a double-bubble case, the number of transition frequencies differs 
from the number of natural frequencies. This result means that the bubbles 
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in the latter case have a characteristic frequency that differs from the natural 
frequency. 

First, we derive a mathematical model to use for the present investigation. This 
model describes the pulsation of two coupled bubbles. Since we are interested 
only in the qualitative natures of the transition frequencies, we use a classical 
theoretical model, called a coupled-oscillator or self-consistent model [10-12], 
although models with greater accuracy have been proposed [4,6,7,9]. 

A gas bubble immersed in a liquid pulsates when a sound wave is applied. 
The sound pressure at the bubble position drives the pulsation. When other 
bubbles (named "bubble 2" ~ "bubble N" , where N is the total number of 
bubbles) exist near the bubble (bubble 1), the sound waves scattered by the 
other bubbles also drive the pulsation of bubble 1. Namely, the driving pressure 
acting on bubble 1, p<n, is expressed as 

N 

Pdi =Pc X + J2PsiJ> (!) 

J'=2 



where p ex and p s ij are the sound pressures of the external sound field and the 
scattered wave emitted by bubble j, respectively, at the position of bubble 
1. When the wavelength of the external sound wave is much larger than the 
bubbles' radii and the distances between the bubbles, p cx can be considered 
spatially uniform. By assuming that the sphericity of the bubbles is thoroughly 
maintained and that the liquid surrounding the bubbles is incompressible, the 
scattered pressure can be estimated with (see, e.g., [15]) 



where p is the density of the liquid, r\j is the distance between the centers of 
bubble 1 and bubble j, Rj is the radius of bubble j, and the over dot denotes 
the time derivative. 

The radial oscillation of a spherical bubble in an incompressible viscous liquid 
is described by the RPNNP equation [16]: 

3 - 2 1 1 

RiR 1 + -R 1 Pwi = — Pdi, (3) 
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where P is the static pressure, a is the surface-tension coefficient at the bubble 
surface, ji is the viscosity of the liquid, k is the polytropic exponent of the gas 
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inside the bubbles, and R\o is the equilibrium radius of bubble 1. Substituting 
Eqs. (1) and (2) into Eq. (3) yields 

3 1 1 ^ d 

RiRi + -Rf -- Pwl = —box + E —-jiitfRj)]- 
l p p ■ 2 T\j at 



When N = 2, this equation is reduced to 



3 1 1 d 

RiRi + -R\ ~-Pwi = ~[Pex + ^j t {RlR2)\, (4) 



where D = ri 2 (= r 2 i). Exchanging subscripts 1 and 2 yields the model 
equation for bubble 2: 

R2R2 + \r\ - -P W 2 = --box + ^(R\Ri)l (5) 
2 p p JJ at 
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By assuming that the time-dependent bubble radii can be represented as R± = 
Rio + ei, R 2 = R20 + e-2, and |e x | < i?i , \e 2 \ < i? 2 o, Eqs. (4) and (5) are 
reduced to the following linear formulae: 

ei + u l0 ei + 5iei = ~^-^e 2 , (6) 

priio tx\QL> 

e 2 + cj 20 e 2 + <) 2 e 2 = D n ei, (7) 



where 



'jo 



1 



pi? 2 



■30 



2a 



3kP + (3k - 1) — 



.70 



for j = 1,2, 



are the partial natural (angular) frequencies of the bubbles, and 

4u 

6 i = ^~ for j = 1,2, 

are the damping coefficients. (In general, the damping coefficients are deter- 
mined by the sum of viscous, radiation, and thermal damping [17]; we however 
use the above setting to simplify the following discussion.) 
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The system of equations (6) and (7) corresponds to that derived by Shima 
when a — 0, 5i — 0, and 5 2 = [2], and to that by Zabolotskaya when a = 
[5]. This kind of system of differential equations, called the coupled-oscillator 
model or self-consistent model, has been used repeatedly to analyze acoustic 
properties of coupled bubbles [2,5,8,10-13], and has been proved to be a useful 
model that can provide a qualitatively accurate result despite its simple form. 

Let us analyze in detail the transition frequencies of the linear coupled system. 
By assuming that the external sound pressure at the bubble position is written 
in the form of p ex = —P a sin cut, a harmonic steady-state solution of the linear 
coupled system is given as 

ei = K\ sm(u)t — 0i), (8) 



where 



K ' = i^^FTi?. ( 9 ) 



= tan 



-1 (Bx 



with 



_ H lF + M 2 G 

_ H 1 G-M 2 F 
Bi ~ F2 + G2 , (11) 

F = L,L 2 - - Ml M 2 , 

G = L X M 2 + L 2 Mi, H 1 = L 2 + ^V, 

L i = (^10 - ^ 2 )> L 2 = (^20 - ^ 2 )> 
Mi = 61U, M 2 = 5 2 uj. 



Here the solution for only bubble 1 is shown. The phase difference <p\ becomes 
% / 2 when 

A x = 0. (12) 
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It should be noted that a case in which both A\ and B\ become zero does not 
exist. From Eqs. (10) and (11), one obtains 



„ 9 „ 9 H 2 + M 2 



F 2 + G 2 \R loP/ 



The numerator of this equation always has a nonzero value, since M 2 > 0; 
this result denies the existence of a case where both A 1 — and B\ — are 
true. Also, it should be noted that F 2 + G 2 appearing in the denominator of 
Eq. (10) always has a nonzero value. When G = 0, F is reduced to 

M 2 2 R10R20 ,4 w w 



This has a nonzero, negative value because M 2 L 2 /Mi > 0, R\qR2qoj 4l / D 2 > 0, 
and MiM 2 > 0. This result means that no case exists where both F = and 
G = are true. As a consequence, Eq. (12) is reduced to 

H X F + M 2 G = 0. (13) 



In the following, we analyze this equation. 

When the damping terms in Eq. (13) are negligible (but exist), one can easily 
obtain the solution for this equation. By assuming that M x ~ and M 2 ps 0, 
one obtains 

H X F ss 0. 



This equation can be rewritten into two independent equations: 

F L X L 2 — — to = (14) 



and 



H 1 = L 2 + ^cu 2 = 0. (15) 



Equation (14) is the same as the theoretical formula given in Refs. [2,5] to de- 
rive the natural frequencies of a double-bubble system; namely, the transition 
frequencies given by this equation correspond to the natural frequencies. This 
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equation predicts the existence of two natural frequencies, 



^10 



> , + < ± VKo - ^ 2 o) 2 + 4^> 2 2 ^F 

^ = 2(1-^) ' (16) 



and is symmetric; namely, it exchanges subscripts 1 and 2 (or 10 and 20) 
in this equation to reproduce the same equation. This means that the two 
bubbles have the same two natural frequencies. One of the natural frequencies 
converges to u 10 and the other to cu 2 o f° r D — > oo, and the higher natural 
frequency rises and the lower one falls as D decreases [2,18]. 

The solution of Eq. (15) given for the first time is 

2 = ^20 (17) 

1 1-R 20 /D 1 ; 



This converges to io\q for D — > oo, and increases as D decreases. In contrast to 
Eq. (14), Eq. (15) is asymmetric (H\ ^ H2); namely, this serves to break the 
symmetry of the natural frequency mentioned above. The transition frequency 
given by Eq. (17) is not the natural frequency because it cannot be given by 
the natural frequency analysis performed in Refs. [2,5]. (Even the other models 
used in, e.g., [6,8] give only two natural frequencies.) 

The results given above show that, when the damping effect is negligible and 
the radii of bubbles are not equal, the bubbles have three asymmetric tran- 
sition frequencies; one of these, for D — > 00, converges to the partial natural 
frequency of a corresponding bubble, while the remaining two converge to that 
of a neighboring bubble. (In the following, the former transition frequency 
is called "fundamental transition frequency" (FTF), and the latter two are 
called "sub transition frequencies" (STFs).) One of the STFs always increases 
as bubbles approach each other. The other STF decreases (increases) and the 
FTF increases (decreases) when the partial natural frequency of the bubble is 
higher (lower) than that of the neighboring bubble. 

Figure 1 shows the transition frequencies for \i m as a function of I = 
D/(R 10 + i?2o)- Other parameters are set to p = 1000 kg/m 3 , P — 1 atm, 
a = 0.0728 N/m, and k = 1.4. Equations (16) and (17) are used in plotting 
those graphs. The radius of bubble 1 is fixed to 10 fim and that of bubble 2 
varies from 5 /im to 10 /im. As discussed above, three transition frequencies 
changing with D appear in each graph, except for the case of R\o = R20, where 
only one decreasing transition frequency appears. (When the two bubbles have 
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the same radii and pulsate in phase, Eq. (6) is reduced to 

+ — je 1+ u 10 e 1 + 5 1 e 1 = -—. 

This equation has only one transition frequency [5,10] of 

which converges to uj\ q for D — > oo and decreases with decreasing D.) It is 
interesting to point out that the larger bubble, whose partial natural frequency 
is lower than that of the smaller bubble, has the highest transition frequency 
among all the bubbles. 

Now we present numerical solutions of Eq. (13) to examine the influences 
of viscosity on transition frequencies. Figures 2—4 show the results obtained 
by using /i = 1.137 x 1CT 3 kg/(m s), which corresponds to the viscosity of 
water at room temperature. From those figures, we can observe that, as the 
viscous effect grows strong, i.e., the bubbles' radii become smaller, the STFs 
vanish gradually from the large-distance region (and sometimes from the small- 
distance region), and only the FTF remains. In the case of R\q = 1 /im, 
the STFs of a larger bubble disappear, and, when the bubbles are of similar 
size to each other, the FTF and the higher STF of a smaller bubble vanish 
in the small-distance region. In the case of R 10 = 0.1 /im, it is difficult to 
distinguish the STFs from the FTF since only a smooth curve, decreasing 
almost monotonically with decreasing /, appears. (The transition frequency of 
the smaller bubble remaining in the small-distance region may be the lower one 
of the STFs.) Those results show that the influence of the mutual interaction 
of the bubbles on the transition frequencies depends on the viscosity of the 
surrounding material, and that this interaction weakens as the viscous effect 
strengthens, i.e., the threshold of the distance for the appearance of the STFs 
is shortened. 

Lastly, we briefly discuss the pulsation amplitudes of the bubbles. Figure 5 
shows otj = (P a / Rjop)Aj and j3j = (P a / Rjop)Bj for Ri = 10 /im, R 2 = 8 
/im, I = 3, and /i = 1.137 x 10~ 3 kg/(m s) as functions of uj/ujiq, where we 
set P a = O.lPo- The present setting is equivalent to that used for one of the 
graphs shown in Fig. 2 except for I. In each graph shown in Fig. 5, we observe 
three reversals of the sign of otj and only two resonance responses, as expected. 
At the transition frequency that does not correspond to the natural frequency 
(the highest one of bubble 1 and the second highest one of bubble 2), no 
resonance response is obtained. This result confirms that the characteristic 
frequency given by Eq. (15) is not the natural frequency. 
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In summary, it was predicted in this Letter that a gas bubble interacting acous- 
tically with a neighboring bubble has three transition frequencies that make 
the phase difference between the bubble's pulsation and an external sound 
be 7r/2, while readymade theories predict only two natural frequencies. This 
present result shows that, in a double-bubble case, the phase reversal of a bub- 
ble (e.g., from in-phase to out-of-phase with an external sound) can take place 
not only at the bubble's natural frequencies but also at some other frequency. 
The present theory for transition frequencies may be useful for understand- 
ing the reversal of the sign of the secondary Bjerknes force, which has so far 
been interpreted using only the natural frequencies [5,19-21]. (This subject is 
discussed in a separate paper of ours [22].) Furthermore, the present results 
might also affect understandings of some related subjects such as acoustic lo- 
calization [23,24] and superresonances [10], previous investigations of which 
have employed systems containing only identical bubbles. The present theory 
for a double-bubble system will be verified [25] by the direct numerical sim- 
ulation technique proposed in Ref. [26], and will be extended to a theory for 
an iV-bubble system [27], where N denotes an arbitrary positive integer. 
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Fig. 1. Transition frequencies of bubbles 1 (uj±) and 2 (0)2) for Rjo ~ 10 (im 
and \i ps 0, normalized by wio- The dashed line denotes Uj/uio = 1. The highest 
transition frequency of bubble 1 and the second highest one of bubble 2 are given 
by Eq. (17). 
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Fig. 2. Transition frequencies for Rjo ~ 10 fim normalized by u>io- The viscous 
effect is taken into account. The dashed line denotes ujj/uj\o = 1. 
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Fig. 4. Same as Fig. 2 except for Rjq ~ 0.1 fj,m. 
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Fig. 5. ctj = (P a /Rjop)Aj and (3j = (P a /Rjop)Bj for i?i = 10 /tm, R2 = 8 /im, and 
/ :!. 
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